Temperature dependence of the charge carrier mobility in gated 
quasi-one-dimensional systems 



CO 

o 
o 

CN 



co 



CZ3 



G 

O 

o 



Lazaros K. Gallos 
Department of Physics, University of Thessaloniki, 54124, Thessaloniki, Greece 

Bijan Movaghar 

SOMS Centre, Department of Chemistry, University of Leeds, Leeds LS2-9JT, United Kingdom 

Laurens D.A. Siebbeles 

Interfaculty Reactor Institute, Delft University of Technology, Mekelweg 15, 2629 JB Delft, The Netherlands 

(Dated: February 2, 2008) 

The many-body Monte Carlo method is used to evaluate the frequency dependent conductivity and 
the average mobility of a system of hopping charges, electronic or ionic on a one-dimensional chain 
or channel of finite length. Two cases are considered: the chain is connected to electrodes and in the 
other case the chain is confined giving zero dc conduction. The concentration of charge is varied using 
a gate electrode. At low temperatures and with the presence of an injection barrier, the mobility 
is an oscillatory function of density. This is due to the phenomenon of charge density pinning. 
Mobility changes occur due to the co-operative pinning and unpinning of the distribution. At high 
temperatures, we find that the electron-electron interaction reduces the mobility monotonically with 
density, but perhaps not as much as one might intuitively expect because the path summation favour 
the "in-phase contributions" to the mobility, i.e. the sequential paths in which the carriers have to 
wait for the one in front to exit and so on. The carrier interactions produce a frequency dependent 
mobility which is of the same order as the change in the dc mobility with density, i.e. it is a 
comparably weak effect. However, when combined with an injection barrier or intrinsic disorder, 
the interactions reduce the free volume and amplify disorder by making it non-local and this can 
explain the too early onset of frequency dependence in the conductivity of some high mobility 
quasi-one-dimensional organic materials. 

PACS numbers: 72.20.-i, 61.30.-v 
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I. INTRODUCTION 

The purpose of this paper is to understand how 
electron-electron interactions affect the ac conductivity 
in the stochastic regime of transport in low dimensional 
materials. Transport in one-dimensional or quasi one- 
dimensional molecular wires formed with self-assembled 
discotic moleculesi*2i2i2i£i§i£, nanotubea&tl at high tem- 
peratures or with disorder, and MBE gate engineered 
confined charged channels, constitutes currently a very 
hot and technically important topic 10 . Other important 
systems include 2- dimensional layers such as those en- 
countered in smectic liquid crystalsii*^, gate charge in- 
jected or doped polymers^, gate addressed FET devices 
based on organic materials^, thin film transistors, and 
same such systems when they exhibit a metal insulator 
transition and superconductors. 

One-carrier conduction phenomena are now reasonably 
well understood, and many workers are in the process of 
designing structures in order to exploit the subtle phe- 
nomena associated with one-carrier resonant tunnelling 
through large molecule»i£*i£'i£. The molecules have en- 
ergy levels structure which can be put in and out of res- 
onance to the energy of the incoming waves. This can be 
done using gates and nano-gates. The changes in these 
energy levels caused by contamination or interaction or 
external fields can then be recorded as a new current- 



voltage response and used in various device applications. 
At the same time it is also well known that small low- 
dimensional structures are more susceptible to electron- 
electron interactions. As long as one is dealing with a 
few steady state carriers which can escape before oth- 
ers arrive, mean-field approaches arc in order. Simple 
minded mean-field approaches will however for obvious 
reasons not work in low-dimensional systems in which 
carriers interfere seriously with each others pathways. 
Interactions are known, for example, to play a role in 
carbon nanotubes at low temperatures, and we believe 
that many ac response phenomena in conjugated poly- 
mers and molecular wires have in the past been wrongly 
attributed to pure disorder effects. In reality one is often 
dealing with a combination of both electron-electron in- 
teractions and disorder and it is not easy to differentiate 
the two experimentally unless one has the right theo- 
retical tools. Electronic interactions also come into play 
with high charging, such as occur for example in electrode 
or gate charged FET organic interfaced and thin film 
transistors. The same applies to highly doped materials. 
Charging energy is particularly important in high-band 
gap, low dielectric constant and low-dimensional molecu- 
lar materials such as chemically or injection doped smec- 
tic and discotic liquid crystalsiaSiMiii. 

Charges can be injected into organic materials chemi- 
cally or by high fields and by light and one can measure 
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FIG. 1: Schematic diagram of the gated model system. A 
molecular wire is first charged by a gate, and the conduction 
takes place along the molecular column between the source 
(S) and the drain (D). 

the frequency dependent conductivity. The frequency de- 
pendent conductivity is in principle one of the easiest 
ways of seeing the effect of interactions. The long range 
coulomb forces produce, naively speaking, disorder-like 
potentials which confine and scatter the carriers. They 
generate in principle a frequency dependence on the mo- 
tion which looks similar to the one produced by true dis- 
order. The problem in practice is, however, that one 
always has some disorder and some electrode polarisa- 
tion effects, and then it is difficult to disentangle these 
processes from each other. 

II. THE COMPUTATIONAL MODEL 

Using many particle Monte-Carlo simulations, we 
study the effect of charge-charge interactions on the 
ac conductivity of finite and electrode addressed one- 
dimensional chains following the work in reference El 
A schematic representation of the system is shown in 
figure n We consider that a molecular columnar wire, 
which does not interact with its neighboring columns, 
is charged via a gated electrode and is attached to two 



metal electrodes, which serve the purpose of source and 
drain. The gate electrode lies at a distance of d=l run 
from the column, and its distance from the edge elec- 
trodes equals Lq=1A nm. We have also considered the 
more realistic situation that the gate electrode is d=30 
nm away from the column. A number N of charges (neg- 
ative or positive, in our case negative) is placed along 
this linear gate, which is parallel to the organic column. 
Charges (which can be either holes or electrons, and in 
our case we consider them to be holes) are created in the 
column and move along its length under the influence of 
an externally applied field, which can be either DC or 
AC. These holes can enter or leave the column through 
the metal electrodes, but asymptotically they establish a 
steady-state distribution in the column. 

In practice, we consider a column of 300 sites, with a 
lattice constant a — 0.35 nm, which corresponds to a 0.1 
micron device. The number N of charges in the gate is 
fixed by the assumed value of the gate field near the plate 
right at the start and remains constant throughout the 
simulation. We assume that the negative voltage at the 
gate causes an equal number of charges (holes) to appear 
in random positions in the column. 

The Coulomb energy V n (t) of a hole located on site 
with index n at distance r n = na from the source elec- 
trode, is the sum of repulsive forces due to the presence of 
other holes and the attractive force due to the negatively 
charged gate: 

V n (t) = V n<Iep {t)+V nMt = — — V +K,att , 

m—1 1 1 

(i) 

where H m (t) is the hole density of site m at time i, and 
assumes the values of 1 or 0, depending on whether this 
site is occupied by a hole at time t. 

For the calculation of the Coulomb attraction energy 
on the holes we use the jellium model. We consider that 
the number of charges in the linear gate are iV, the length 
of the column is L, the gate electrode starts at a distance 
Lq from the source and drain electrodes and its distance 
from the column is d. Under these conditions, the attrac- 
tion energy for a hole located on site r n is equal to 
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During each step, the transition probabilities for all the charges in the column are computed, so that a charge in 



3 



site n has a hopping transition rate W ni7l ±i towards its 
neighbors 



W n , n ±i = 



cxp (- AE, ^ l±1 ) iorAE, hn±1 > 
forA£L n± i < 



(3) 

The energy difference between the two sites is computed 
by the formula 



AE n , n ±x = K,n±i - K, T eFa , 



(4) 



where the symbol V n , n ±i represents the Coulombic en- 
ergy on site nil when we displace the hole located on 
site n to the site n ± 1. The jump frequency vq is fixed 
to 10 12 Hz. The external electric field F can either be 
constant Fq —20 kV/cm, or oscillating with a period w 
so that F(t) = Focos(uit). When a neighboring site is 
occupied the transition rate towards that site is 0, i.e. 
we use excluded volume rules. 

For the first and last site of the column the quantity 
AE is modified, since now there is only one neighbour 
to move to, and injection/absorption to the electrodes is 
possible. There exists an energy barrier between the va- 
lence band of the organic material and the Fermi level of 
the metal electrode, whose height we denote by E b . The 
energy difference for a hole injection from the electrode 
to the first column site is 



AE 1 ^ =E b + V 1 - eFa , 



(5) 



and for hole injection from the other electrode to the 
opposite edge is 



AE™± T =E b + V N + eFa . 



(6) 



In the same way, the energy difference for absorption of 
a hole at n = 1 to the adjacent electrode is 



AEf^ =-E b -Vi + eFa 



(7) 



while the same quantity for a hole leaving the column at 
the other side is given by 



A££^ L = -E b -V N - eFa 



(8) 



The case of E b = oo corresponds to blocking barriers, 
where the number of holes in the system is fixed for the 
entire duration of the simulation and no hole can enter 
or leave the column. 

Once the transition probabilities W have been com- 
puted for all holes in the system, the time St required for 
each hop is drawn by an exponential distribution, 



St = 



In J? 



W n ,n-1 + W n ,n+l 



(9) 



with R being a uniformly distributed random variable, 
assuming values between and 1. The shortest time 
<^min of this ensemble (including the time needed for in- 
jection/absorption) is chosen and this jump takes place. 



The total time advances by the same amount St lri \ n , and 
the direction of the jump is decided with a probability 



P, 



n,n±l 



n,n±l 



w, 



n,n+l 



(10) 



When a hole moves from a site to another, the work done 
by the electric field during the time interval St m [ n is cal- 
culated by 8 A — ±eaF(t). The total work done over 
a time interval is added to the total work AA for that 
interval. 

The above procedure is repeated many times with the 
holes moving around, leaving and entering the column. 
The time that a given realization of the system lasts de- 
pends on the convergence of the results, which in turn 
depends on the specific parameters of the problem. Thus, 
typical times simulated are in the range of 0.01-0.1 ms. In 
order to eliminate any initial distribution bias, we repeat 
the whole algorithm for a number of different realizations 
(typically of the order of 10-100), with different random 
initial configurations. During all these realizations we 
monitor the average of the holes mobility /i in the sys- 
tem with time, until this average is well converged. In 
almost all cases, the error bars for the mobility estima- 
tion were less than 10%. 

The number of holes Nh in the column is a number 
which varies with time and converges asymptotically to 
an average value, depending on the exact simulation con- 
ditions. When we use blocking electrodes, of course, Nh 
is always equal to N, the number of negative charges in 
the gate. 

The average hole mobility fi is calculated via the work 
AA done on all holes during an oscillation period At = 
2it/lu, and is given byi£ 



2AA 



eFJAtNh 



(11) 



III. RESULTS 

We have simulated the gated system described above 
for both DC and AC externally applied field. When we 
considered an AC field, we also studied the following two 
cases: a) include a finite barrier height between the metal 
electrodes and the column, and b) assume an infinite bar- 
rier height, which corresponds to blocking electrodes. 



A. DC Mobility 

In figure |21 we present the average mobility for the 
gated one-dimensional column with 300 sites in a DC field 
F=20 kV/cm, as a function of the number of charges N 
in the gate. The injection energy barrier height is equal 
to E b =0.2l eV. Neglecting all carrier-carrier interactions, 
including those from the gate charge, the mobility would 
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FIG. 2: Average carrier mobility as a function of the number 
of charges N in the gate for a DC field of 20 kV/cm. Different 
temperatures are plotted(top to bottom): 450 K, 300 K, 150 
K, and 100 K. Lines are used as guides to the eye. 

scale as 

fj, = /i exp (~^) ■ ( 12 ) 

Different curves correspond to runs at different tempera- 
tures. Figure 2 shows that the temperature dependence 
of the mobility is strongly dependent on the value of N, 
the charge density. The increase in mobility with tem- 
perature is due to two factors, a) The carriers find it 
easier to come into the chain, and as the temperature 
increases because of the barrier factor (i.e. as in equa- 
tion b) the charges can also hop more easily inside 
the chains against the self-field produced by the other 
carriers, which means that in equation 1121 the factor fiQ 
should be replaced by no(T, N). The calculated low tem- 
perature mobility is very small at T=50 K (not shown). 
In this limit the charges are practically immobile and can- 
not overcome the charge-charge interactions which keep 
them vibrating around their "equilibrium" positions. 

An extremely interesting phenomenon is observed 
when we look at the temperature dependence of the cal- 
culated mobility at different densities. Figure [5] shows 
that there are densities at which the "activation energy" 
is practically zero. This means that for these particu- 
lar values of N, a charge distribution is established in 
the given field geometry, where the injection barrier is 
compensated by the internal space charge and the gate 
field. One can also look at it as the variation of the calcu- 
lated mobility with the number of charges N. When the 
temperature is constant, the average mobility decreases 
with increasing N. This is expected because the carriers 
have less space available to move in. This observation 
is actually true for room temperatures and above, but 
figure |21 shows that this rule is clearly violated at lower 
temperatures, and for high concentrations. For low tem- 
peratures, the mobility decreases with charge concentra- 
tion, but when the concentration exceeds a certain value 



(around 7V=30, corresponding to a concentration of 0.1) 
the mobility increases with increasing concentration un- 
til A/=100. The mobility decreases again, but exhibits 
another peak at N=150. For concentrations higher than 
0.5 the average mobility drops extremely rapidly, and 
the charges are practically immobile. In other words, in 
the mobility versus A~ plot, we can observe three max- 
ima (located roughly at N=5, N=100, and N=150) and 
three minima (for N=30, N=120, and AT=180, which 
is the highest concentration we used). One can under- 
stand this behavior if one notes that the combination of 
the intrinsic (work function) barrier, the charging energy 
coming from other carriers, and the charging energy from 
the gate within a particular gate geometry, together pro- 
duce a density dependent pinning barrier which modu- 
lates the absorption and injection efficiency. One can see 
that when absorption/injection efficiency is high, the mo- 
bility can be high, and this is caused by the forward drift 
of all carriers. As soon as the density exceeds the criti- 
cal number N c the injection is reduced and the incoming 
carrier has to wait until one has exited . To calculate N c 
we note that the Coulomb barrier at entry on column is 
defined by 

N 2 

^ 47ree o n(a(A0) 

where (a(iV)) is the average distance between two car- 
riers at total number N. With e ~ 3, Eqb exceeds kT 
when the number of carriers A" on the chain is larger 
than N c ~ 40 — 50 (the value increases slightly with the 
temperature). 

When the Coulomb repulsion energy at the entry site 
starts to get much bigger than kT the carriers block each 
other along the column. They can only move sequen- 
tially, "in a phase" forced on them by the internal fields. 
The total energy barrier seen by the carrier on the first 
site in the column is given by 

^=E b + J2^--E mts , (14) 

where E gatc is the "attractive" energy of the gate counte- 
rions and depends on the distance and gate geometry, and 
the sum runs over the position of the carriers. When we 
increase the distance to the gate counter charge distribu- 
tion from 1 nm to a more realistic value for conventional 
gate devices of d = 30 nm (see figure |5|l 1 4 i 19 i 20 , the mo- 
bility versus density behaves quite differently as shown 
in figure |21 The internal charge density as a function 
of gate voltage is also shown in figure [3J Note that we 
now have fewer charges for the same gate charge because 
the coulomb repulsion is relatively stronger and prevents 
charges from coming into the column. Note also that 
the mobility versus A~ oscillations observed in figure [5] no 
longer appear in the same place. 

For higher temperatures, the picture is much easier to 
understand because the charges have enough energy to 
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FIG. 3: Comparison between systems with d=lnm (filled 
symbols) and d=30 nm (open symbols) of a) the average mo- 
bility and b) the average number of charges in the column, 
for T=100 K (circles) and T=300 K (squares), under the in- 
fluence of a DC field of 20 kV /cm. 



overcome the local attractive energy barriers and there 
is no pinning effect. When T=450 K and N—5, for ex- 
ample, it is much easier for holes to enter the system, 
so that there are always 8 to 9 charges in the column 
which undergo a forward drift motion, as opposed to 1 
or 2 charges in the column for T=100 K. 

For the injection barrier Eb value used, the average 
number of holes in the column during the simulation is 
close to the value expected from the gate charge, except 
of course when we use d=30 nm. Here, the Coulomb re- 
pulsion with weak screening keeps the charges out of the 
column (see figure 01. For d=l nm, the gate to charge 
density matching is almost one to one for all tempera- 
tures, and at the lower end of concentrations. At higher 
concentrations, we observe a reduction in the average 
number of holes in the column of the order of 10-15%, 
e.g. for iV=100 we have around 90 holes. This is not sur- 
prising and is an effect of the "not metallically screened 
coulomb repulsion" in the chaini2i2£. 
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FIG. 4: Temperature dependence of the average AC mobility 
for different initial charge concentrations and field frequen- 
cies: a) 1 6 Hz, b) I0 9 Hz, and c) I0 10 Hz. The energy barrier 
height of the electrode is equal to 0.21 eV. We observe differ- 
ent patterns of behavior as a function of the concentration as 
we vary the frequency. The gate charges are as follows: N=3 
(□), iV=9 (<), AT=18 (•), 7V=36 (0), and iV=60 (x). 
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B. AC Mobility 

1. Finite electrode barrier height 

For the relatively low frequencies of 10 6 Hz (figure^), 
and for intermediate densities, the mobility increases al- 
most one order of magnitude as we increase the tem- 
perature (the injection barrier is limiting the mobility). 
But, at the low concentration end, when carrier interac- 
tions are small the temperature dependence is weak. At 
very low frequencies, one should of course expect similar 
results to the ones observed in the limit of a DC field. 
So, the temperature independence at low concentrations 
implies that at 10 6 Hz one is already looking at bulk mo- 
tion. The barrier no longer affects the "single-particle" 
mobility. It is important to note that apparently all the 
temperature dependence in figure 2t is already due to 
many-body effects. When the frequency of the electric 
field is increased further, as shown in figure 2Jd, to the 
frequency value of 10 9 Hz for example, the change with 
temperature is now weaker even for higher densities. On 
this time scale, the interactions are getting less important 
as one might expect. Indeed, for very high frequencies, 
e.g. 10 10 Hz as shown in figure^, the mobility actually 
decreases with increasing temperature, and this behavior 
is monotonic with N . The change in mobility is however 
less than half an order of magnitude. This behavior is 
unusual for hopping systems^. We believe that this is 
because the effective disorder is increasing with temper- 
ature. The disorder increases because the number of car- 
riers in the column is increasing and this reduces the free 
volume, as shown in figure [21 

For a given temperature, the calculated frequency de- 
pendence is shown in figure |S] Here we have plotted 3 
curves: a) T=193 K, b) T=343 K, and c) T=493 K. Con- 
sider first the low temperature curve. Here the average 
mobility at small concentrations is practically constant 
for all frequencies. This means that the injection barrier 
has very little effect on the frequency behavior at this 
temperature. Let us try to understand why. In the spirit 
of a simple resistance sum2£ for a one particle random 
walker and rigorous in the low frequency limit, we have 
for a column of 300 sites a transport admittance Y: 



1 



299 



kT \icu + W h 



iuj - 



(15) 



Wbulk , 

The injection rate W ln ^ is activated as in equation 1121 
and the bulk rate Wbuik is a constant for a single parti- 
cle. From equation 1151 one would expect the frequency 
dependence to come in when the frequency is of the or- 
der of the injection rate, that is roughly 10 8 Hz, and to 
saturate at the bulk value of 10 12 . This, however, does 
not happen in figure El At N = 3 the mobility stays 
roughly at saturation value all the way. The reason for 
this behavior is as shown in figure [2] The particles in 
the simulation are moving in a potential made up of the 
gate field, the electron-electron interaction and the bar- 
rier, i.e. as given by equation 1141 The barrier at a given 
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FIG. 5: Variation of the average AC mobility with the field 
frequency. The energy barrier height of the electrode is equal 
to 0.21 eV. Three different temperatures are presented: a) 
T=193 K, b) T=343 K, and c) T=493 K. Symbols correspond 
to: iV=3 (□), JV=18 (•), iV=36 (A), iV=60 (x). 



concentration may indeed be much lower than the injec- 
tion barrier, and that is why there is no, or only a weak 
frequency dependence for this particular concentration. 
Looking at figure we see that at higher concentrations 
the mobility becomes strongly frequency dependent (see 
the corresponding dc limit in figure EJ), and rapidly in- 
creases with frequency at around 10 9 Hz. This is now an 
interaction effect as in equation 1141 and it looks as if in 
the spirit of one body theory, interactions are causing an 
effective disorder. The interactions appear to have little 
effect on time scales shorter than 10~ 9 s. As we raise the 
temperature further to 343 K, the frequency dependence 
is now weak for all N. Above 343 K, there is in effect no 
strong inhomogeneity in the hop time distribution. This 
means that the variation in hopping activation energies is 
comparable to kT . Finally, in the extreme high T limit, 
when T=493 K the mobility actually decreases with in- 
creasing frequency. This is most unusual for hopping 
transport and is in this model due to the extra charging 
that high temperature causes, as shown in figure [3 At 
high densities, a real time analysis shows that the charges 
move like a collective "line of billiard balls" . A carrier 
can only come in when one has escaped, and the mobility 
decreases with N at high temperatures. In figure we 
present the steady state number of holes in the column 
as a function of temperature for different gate charges. 
This number is an average over all different frequencies 
used. It is evident that there is a considerable net charg- 
ing effect in the system. The charging is more prominent 
for smaller N. For example, when 7V=3 and T=593 K 
the average number of holes in the system is around 7, 
while for T=193 K it is 2 (the absolute value of these 
numbers depends of course on the value of the interface 
energy barrier Eb, but the basic trends remain the same). 
This charging effect is responsible for giving the lower 
value of the mobility observed at higher temperatures 
(see for example figure [5| • Normally, in hopping trans- 
port we expect the mobility to go up with temperature. 
Here however, the steady state density of holes, which 
also influences the magnitude of the hole-hole scattering, 
is temperature-dependent. There are competing influ- 
ences: whereas temperature makes hopping easier, the 
number of holes also goes up, and this increases there- 
fore the hole-hole interactions, or what is also the degree 
of disorder in the system. 



2. Blocking electrodes 

We also studied the case where there is an infinitely 
high energy barrier among the electrodes and the molec- 
ular wire. This corresponds to blocking electrodes, which 
means that the system is "closed" in the sense that no 
carrier can (subsequent to charging) enter or leave the 
column. In figure we can see that at high frequencies 
the average mobility is constant and almost the same for 
different concentrations. This is because in this limit we 
are looking at bulk hopping. As we lower the frequency, 
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FIG. 6: Average over all AC frequencies of the number of 
holes along the column, as a function of temperature. The 
energy barrier height of the electrode is equal to 0.21 eV. Top 
to bottom: A?=60, A=36, N = 18, and N = 3. Note the 
chareine effect at hieher temDeratures and lower concentra- 
t 
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FIG. 7: Average mobility as a function of frequency for dif- 
ferent charge concentrations (top to bottom: N=3, N=18, 
N—36 ) with blocking electrodes and 300 sites at T=293 K. 



the mobility drops steadily to very small values. The 
carrier displacement is at low frequencies limited by the 
boundary and the decrease in mobility is expected in a 
confined system. The effect is also concentration depen- 
dent. Going up in frequency, we observe that the higher 
the hole densities, the later the saturation frequency. For 
the same frequency, the difference in the mobility for the 
low and high concentrations can be up to 2 orders of 
magnitude. This is true for concentrations in the range 
from 0.01 to 0.12 charges per site. 

The corresponding temperature/frequency dependence 
is shown in figure [SJ The curves seem to lie practically all 
on the same points except in the limit of high frequencies. 
Here, the average mobility is temperature dependent and 
the interesting point is that the mobility decreases with 
temperature. This is a many-body effect and says that 
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FIG. 8: Average mobility as a function of frequency for dif- 
ferent temperatures with blocking electrodes and 300 sites, 
where N=36. Symbols correspond to the following tempera- 
tures: T=193 K (•), T=293 K (□), T=493 K (x), and T=593 
K(A). 

temperature is enhancing the coulombic disorder by al- 
lowing configurations which provide a wider distribution 
of activation energies. 

In order to understand this behavior more accurately, 
we monitored the motion of the charges in the column 
as a function of time. The main conclusion is that the 
charges follow the periodic pattern imposed by the exter- 
nal alternating field. When the field frequency is high, 
the charges do not have enough time to come closer to 
each other or to the electrodes, so they fluctuate rapidly 
around their initial positions following the field, and this 
gives rise to the given calculated mobility value shown in 
figures [7] and |21 As the temperature increases, the ther- 
mal energy of the charges allows them to move more eas- 
ily against the direction of the fields. This increases the 
electron-electron disorder energy, and in this way tem- 
perature now reduces the overall mobility! This picture 
explains the decrease of the mobility with temperature 
at a given value of (high) frequency, as shown in figure |S] 

For a given high frequency range, this picture remains 
basically the same. Suppose that 3 charges are located 
in the column at some equilibrium positions. For a fre- 
quency 10 12 Hz, for example, they will make very short 
moves towards one or the other side. When the fre- 
quency is switched to 10 10 Hz, the charges will again 
follow the field, but despite the larger distances traveled 
in any one direction, the charge closest to the electrode 
will still not have enough time to seriously be blocked by 
the confining electrodes. Similarly, the other charges will 
not have enough time to feel the changes in each others 
positions, and they will continue to oscillate, essentially 
undisturbed, thus giving rise to the same mobility as for 
higher frequencies. This explains the constant value of 
the n(N) mobility versus frequency at high frequencies. 

When the frequency is small so that the displacement 
can reach the electrode boundary, the charges are pushed 



towards one side of the column. They come very close 
to each other and thus the mobility decreases with N 
(only the one furthest from the edge can move easily, but 
when the charge tries to escape, the external field pushes 
it back). As the frequency gets lower and lower, the effect 
becomes more pronounced because now the biasing field 
remains at the same magnitude and in the same direction 
for a longer time l/u>. This is the reason why the mobility 
falls rapidly and monotonically with lower frequencies. In 
fact, when using blocking electrodes, the DC mobility is 
rigorously as it should be, with the (+) charges 'piling' 
up near the (-) electrode. 



IV. CONCLUSIONS 

We have continued the study of charge transport in in- 
teracting quasi-one- dimensional gated systems with em- 
phasis on the temperature dependence of the mobility. 
We have seen that the temperature changes the mobility 
in the presence of a DC field but that the effect is very 
much dependent on charge concentration. The temper- 
ature dependence is partially due to the assumed injec- 
tion barrier and partially due to the fact that carriers can 
"push their way in" like a line of "billiard balls". Both 
factors work together to produce the final result shown 
in figure Temperature can cause the net steady state 
charge to vary and this then also affects the charge-charge 
interaction efficiency (figure |2J). 

The many body interactions do indeed cause a fre- 
quency dependence of the mobility. However, using the 
same parameters as for doped discotic HAT6^, we find 
that the frequency dependence comes in at too high a 
frequenc y. It seems that the experimental data in refer- 
ences |l||2||3j^|23j can only be understood if we also include 
disorder. 

When we have blocking electrodes, the interactions in- 
fluence the density and the temperature dependent ac 
mobility. The charges increase the confinement and the 
mobility decreases with concentration and with temper- 
ature in the way one would expect for finite ordered 
segments. Finally, we now need to extend our theo- 
retical tools to handle 2-dimensional systems, so that 
we can model smectic liquid crystals, mixes of smectics, 
FET's, thin-film transistors (TFT), and also the trans- 
port through large molecules such as nanotubes and nan- 
otubes filled with fullerenes (pea pods) 2 ^. So far the work 
has focused on quantum transport. There are and will 
be, we believe, many new and exciting experimental re- 
sults coming up in this field in the stochastic and diffusion 
limit. 
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